Delocalization induced by low-frequency driving in disordered superlattices 
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We study the localization properties of disordered semiconductor superlattices driven by ac-fields. 
The localization length of the electrons in the superlattice increases when the frequency of the driving 
field is smaller than the miniband width. We show that there is an optimal value of the amplitude 
of the driving field for which the localization length of the system is maximal. This maximum 
localization length increases with the inverse of the driving frequency. 
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Real materials always contain a certain degree of disor- 
der since the atomic structure is never perfectly regular. 
In fact, many physical properties are either influenced or 
even mainly determined by this randomness. The under- 
standing of the effects of disorder in the physical prop- 
erties of a material is therefore of great practical impor- 
tance and has played a central role in condensed matter 
physics in the last half-century 0. 

One of the most simple disordered systems is the mo- 
tion of a particle in a one-dimensional random potential. 
Realizations of this system can be studied experimen- 
tally, e.g. in semiconductor superlattices (SL) made by 
varying the alloy composition along one of the dimensions 
of a compound semiconductor such as AlxGai_xAs Q. 
With the technique of molecular-beam epitaxy the de- 
gree of control of the disorder in such systems can be very 
high. If the SL is long enough its eigenstates are exponen- 
tially localized due to disorder. Electrons in these local- 
ized states are spatially confined and their only contribu- 
tion to transport is through thermally activated hopping. 
The main quantity of interest in this case is the local- 
ization length A of the electron wave-functions, which is 
controlled by the ratio between the bandwidth A and 
the strength of the disorder W. A disordered system of 
length L > X will behave as an insulator while a system 
with length L < X will behave as a conductor 

The fascinating effects of radiation on the transport 
properties of SLs have been intensively studied both the- 
oretically and experimentally in the last three decades 
0. One of the most interesting effects, first predicted 
by Dunlap and Kenkre, is dynamic localization (as de- 
fined in 5]), in which the electron wave- function can be 
strongly localized in the presence of an AC electric field. 
This effect was shown to be associated with miniband 
collapse which occurs at the zeros of a Bessel function 
that depends on the field amplitude. At these points the 
width of the miniband goes to zero, the group velocity 
of a wave-packet vanishes and the electron becomes ef- 
fectively localized. Experimentally, this effect was first 
observed as a suppression of current at some amplitudes 
of the AC electric field jTl . 

When a time-periodic driving is applied to a disor- 
der SL, a new possibility emerges, i.e. controlling the 
localization of the system by changing the amplitude 



and/or frequency of the driving field. Holthaus et al. [g 
showed that for high frequencies, the Floquet states be- 
come localized (in contrast with dynamical localization) 
and their localization (participation ratio) depends on 
the ratio A^ff/W, where Ag// = AJo{eVd/Huj). Here, 
V and Lu are the amplitude and frequency of the driving 
field and d is the spatial period of the lattice. 

In this letter we calculate, for the first time, the lo- 
calization length of a driven disordered SL and we focus 
on the previously unexplored low frequency regime. For 
this purpose we use a Floquet-Green function formal- 
ism which makes use of matrix continued fractions '9|. 
This formalism allows us to generalize the definition of A 
for time-periodic systems and to calculate its value. We 
found that there are two distinct and clear-cut regimes in 
this system: High-frequency regime for hoj > (A -I- W)/2 
and low-frequency regime for Hcu < {A -\- W)/2. In the 
first one, we found that A is a function of Aef / /W, in per- 
fect agreement with previous works In this regime, 
the driving field always contributes to localize even more 
the wave- function (as compared to the non-driven case). 
In contrast, in the low frequency regime, we find that 
the driving can significantly delocalize the electrons. In 
addition to this new result, we provide an intuitive ex- 
planation of why this should be so: For low frequency, 
each additional Floquet channel created by the driving 
provides the electron with new paths which will differ in 
their degree of localization, with some having smaller and 
others having greater localization length as compared to 
the non-driven case. Since A is determined only by the 
path with the greatest localization length, after the en- 
semble average has been performed, additional propa- 
gation channels should always contribute to increase A. 
The same situation does not occur in the high-frequency 
case because new paths introduced by the absorbtion or 
emission of one or more photons, always have localization 
lengths smaller than in the non-driven case. 

We will model the one-dimensional disordered SL in 
the presence of an AC field by a single-band Ander- 
son Hamiltonian with diagonal disorder [lOj plus a time- 
periodic potential. 
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The on-site energies Cj are distributed uniformly from 
-W/2 to W/2. The driving potential is due to an AC 
electric field of amplitude 2V. (The factor of 2 is for 
convenience.) In a SL this term can represent THz ra- 
diation linearly polarized in the growth direction of the 
lattice. The quantity A is equal to the bandwidth of the 
non-driven system without disorder. We will use units in 
which h = 1 and we set A = 4. 

For a system that obeys discrete time-translational 
symmetry with period T, there exists a complete set 
of solutions to the Schrodinger equation of the form 
|*"'P(t)) = exp{-iea,,pt/h)\(t>"''P{t)), where |</)"'P(0) = 
\(j)'^'P{t + T)). These are called the Floquet states of the 
system and the periodic functions \(j)"'P{t)) obey an eigen- 
value equation similar to the static Schrodinger equation, 



H{t) - ih 



dt 



ir'''(<))=e«,p|r'^W), 



(2) 



with H{t) the time-periodic Hamiltonian of the system. 
The eigenvalues of this equation can be written as e^.p = 
ta +pfiuj, for < ea < fiLO and p an integer. A Floquet- 
Green operator corresponding to Eq. |21 can be defined 
and its Fourier components satisfy 
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FIG. 1: (Color online) (InGiL(O)) as a function of the length 
of the system L for different values of disorder W and fre- 
quency uj. High-frequency results, uj — 20, show regular os- 
cillations on top of the exponential decay, according to the 
equation G'il(O) = AJo(2VL/huj) exp''^' ^ . Low frequency 
curves do not show these oscillations. The value of A is the 
negative inverse slope of the linear fits. 



the Hilbert space of the system. The Floquet compo- 
nents of the Green operator for this Hamiltonian satisfy. 



4>t+p){(t''v 
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where |</'p°) are the Fourier components of the Floquet 



eigenf unctions 10*^ 



-ipujt 



1 ) . The transport 



properties of driven systems have been formulated in 
terms of this Floquet-Green operator which plays a simi- 
lar role to the Green operator in the Landauer formulism 
for conduction prj. We generahze the usual definition of 
the localization length in terms of Green functions and 
define the localization length of the kth Floquet compo- 
nent as 
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The different quantities g['2 (E) are associated with the 
probability of a process where an electron starts with 
an energy E at site 1 and ends at site L with energy 
E+khu. Each one of these processes in principle can have 
a different localization length associated to it. However, 
assuming that the asymptotic behavior of the Floquet 
wave function is exponentially decreasing, it is easy to 
see that the dominant term in the sum over p in Eq. (O 
will always be the same, independently of the value of k. 
Therefore the functions A^*^^ (E) are identical. From now 
on we will only compute X{E) = \^^\E). 

For the calculation of G^^\E) we use a method de- 
veloped by one of the authors (see Q for details). We 
want to calculate the Floquet-Green operator for a peri- 
odic Hamiltonian of the form H{t) = Hq + 2cos{ujt)V, 
where Hq and V are any time-independent operators in 



(E + khuj - Ho)G^''^ - ^(GC^+i) + G^^-i)) = 4,0. (5) 

These equations can be solved using matrix continued 
fractions. For the case fc = 0, one gets 

G^''\E) = iE-Ho-V,fAE))-\ (6) 



where 



with 



Veff^V+^{E) + V~^f{E), (7) 



K//(i?) = V- 
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E±lhLu-Ho-V- 



E±2hLo- Hq - V—V 



-V 



(8) 



The convergence of equation ||SJ) is system specific. For 
our Hamiltonian, Eq.(l), the number of bands necessary 
to ensure convergence increases linearly with VL/uj. The 
numerical performance of our method is determined by 
the speed in the calculation of an i x i-matrix inverse 
for each Floquet sideband. For our disordered system we 
know that the localization length in the middle of the 

band [E = 0) for = behaves as Ao « 105 

|T^ . In this case it is also known that A as a function 
of the energy E follows a parabolic law with a maximum 
at = 0. For \E\ > (A + W)/2 the locahzation length 
decreases rapidly with \E\. In this work, we will always 
take £: = 0. 

In Fig^we show some examples of the results obtained 
for the ensemble average of In Gil (0) as a function of the 
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FIG. 2: High frequency results. Results for different values 
of disorder W with uj — 20.0. The corresponding results for 
a non-driven system with renormalized bandwidth are shown 
with lines (obtained from numerical data). 



length of the system L. (We have dropped the superscript 
in the Floquet-Green operator.) As expected, a straight 
hne fits the data very weU. The wave-functions decay 
exponentially with the distance. The negative slope of 
these curves corresponds to the inverse of the localiza- 
tion length when L ^ X. It is important to point out 
that the deviations from this exponential behavior are not 
due to a deficient ensemble average. For high-frequencies 
they show a regular pattern, a manifestation of dynamic 
localization which shuts off the tunnelling between the 
sites L — 1 and L at the zeros of Jo{2VL/uj). In fact, 
for high-frequency and L ^ X one can show that the 
Floquct-Grccn function can be expressed analytically as 
Gil(0) = AJo{2VL/u})exp{~L/X). As shown in Fig. [H 
for the second and last sets of data (from top to bot- 
tom), the thin (red online) continuous line representing 
this analytical result fits the numerical data very well. 
For low-frequencies, the Green function also decays ex- 
ponentially. However, the deviations from this behav- 
ior cannot be expressed analytically and seem to occur 
equally on both sides of the straight-line fit. 

We now show results for the high frequency regime, 
which has been discussed in the literature although 
to the best of our knowledge, no calculation of A has been 
reported. The high frequency limit can be characterized 
as the regime in which the absorption or emission of any 
number of photons would leave the particle with an en- 
ergy outside of the region where the eigenenergies of the 
non-driven system concentrate. This region is well known 
to have a width A -t- W, and therefore, at the center of 
the band this condition is satisfied when uj > (A + W) /2. 
In the high-frequency regime, the results for A are ob- 
tainable from the function Xa{A/W) but with a renor- 
malized hopping term A AJ(){2V/uj). Fig. |21 shows, 
for oj — 20.0 and for several different values of disorder, 
that the numerical data is in excellent agreement with 
Holthaus' result. The minima of A correspond to the 
zeros of the Bessel function Jq{2V/uj). 
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FIG. 3: (Color online) A as a function of V/co for W ^ W 
(top panel) and W — 5 (bottom panel) and different values 
of UJ. 



In FigOlwe show the localization length A as a function 
of V for low and high frequencies and for two values of 
disorder, W = 10.0 (top) and W = 5.0 (bottom). Here, 
one can see that there is a fundamental difference be- 
tween those two regimes: For w < (A -I- W)/2 and V 
small, A increases with V, whereas for lu > {A + W)/2 
it decreases. This delocalizatioir of the wave function 
due to low frequency driving is the main result of this 
letter. For all the frequencies that we explored in the 
low frequency regime, as V was increased A initially in- 
creased, then reached a maximum, and finally decreased 
(with some oscillations). This behavior can intuitively 
be understood if one assumes that the driving provides 
the electron with new propagating channels n, each one 
with a different localization length A^. According to this, 
one expects A = max{Xi, X2, ..A„}. Given the linear de- 
pendence of the total number of Floquet sidebands with 
V, it is expected that n (x V, and therefore A should 
initially increase with V. This is valid until we reach 
the maximum number of effective channels that can be 
supported in the region of the spectrum where the local- 
ization lengths are comparable (size w A -I- W). Here, 
A reaches a maximum value Xmax- Beyond this point, 
A decreases with V (with oscillations), as the weight of 
the Floquet eigenstates oscillates and shifts to energies 
outside the band. 

In Fig. 4 we plot the data for {Xmax — Ao)/Ao vs. l/oj. 
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FIG. 4: Behavior of \max as a function of l/u for different 
values of disorder W . The numerical data has been fitted 
using a function \{nmax,(^) derived from a simple statistical 
model where the maximum number of effective channels is 
inversely proportional to the frequency, rimax = cju and a 
is the standard deviation of each channel (with exponential 
probability function). 



This figure suggests that low-frequency driving is more 
effective in delocalizing a weakly disordered system than 
a strongly disordered one. The continuous curves cor- 
respond to a fit to a very simple mathematical model 
where the maximum number of propagating channels is 
inversely proportional to the frequency Umax = cjio. We 
also assume that each one these channels has an expo- 
nential distribution function with standard deviation u 
(a Gaussian distribution gives similar results). We can 
see that this simple model fits the data very well. De- 
spite of the fact that some of the assumptions in this 
simple model are not formally justified, we believe that 
its agreement with many of the features in the numer- 
ical data support our intuitive explanation for the de- 
localization of the wave- function at low frequencies. In 
our system, \max seems to increases monotonically with 
the inverse of w. However, in a real SL a limiting time- 
scale might appear, i.e. the energy relaxation time r. 
If LOT < 1 the delocalization induced by the driving field 
could be limited. For a typical relaxation time r = 10~^^ 



s and a typical miniband width A = 3.5 meV, illumi- 
nation with THz radiation of frequency lu = A gives 
u!T Rs 5 > 1. This implies that a transition between the 
two frequency regimes should be possible to implement 
with THz-irradiated SLs. For a disordered SL of length 
L w Aq a change of several orders of magnitude in the dy- 
namic conductivity should be observed when going from 
frequencies higher than the bandwidth to smaller fre- 
quencies. Also, the behavior of the transport properties 
with temperature should be very different, from a hop- 
ping regime with some activation energy which is favored 
by higher temperatures to a normal regime where trans- 
port processes should be hindered by phonons and other 
temperature related effects. Other systems can also be 
used to observe these delocalizing effects, e.g. atom opti- 
cal traps, which have been proposed as a testing ground 
of Anderson localization. For this system, the lattice 
can be implemented by far detuned counter-propagating 
laser beams, the AC driving can be obtained using a peri- 
odic phase-shift between the beams and the random 
potential can be obtained with a superimposed random 
speckle pattern j3- 

In this work we have calculated the localization length 
for one-dimensional SLs in a homogeneous AC electric 
field. We have found two very different regimes accord- 
ing to whether the frequency of the driving is smaller 
or bigger than the miniband width of the disordered 
SL. For low frequency driving, the localization length 
increases for moderate values of the driving amplitude, 
while for high frequency the localization length decreases 
due to the Bessel function renormalization of the mini- 
band width. We have shown that low-frequency driv- 
ing can increase the localization length by providing new 
propagation channels. In this sense, the effect of such a 
field can be compared to an increase in the dimensionality 
of the system, which suggests that, in disordered systems 
where a metal-insulator transition is expected when the 
dimension of the system increases, low-frequency driving 
could have very dramatic delocalizing effects. 
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